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ABSTRACT 

We have carried out a combined Hubble Space Telescope (HST/GHRS) and 
Far Ultraviolet Spectroscopic Explorer (FUSE) analysis of the prototype dwarf 
nova SS Cygni during quiescence. The FUSE and HST spectra were obtained 
at comparable times after outburst and have matching flux levels where the two 
spectra overlap. In our synthetic spectral analysis, we have used SS Cygni's ac- 
curate HST FGS paral lax giving d = 166 pc, a newly determined mass for the 



accreting white dwarf ( jBitner et al. 1 120071 ) of M^d = O.8IM0 (lower than the 



previous, widely used 1.2Mq) and the reddening (Ep ^y) values 0.04 ( IVerbunt 



1987l : lLaDouslll99lh & 0.07 (IBruch fc Engel I IiQqJ derived from the 2175 A 



^Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space 
Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, 
Inc. under NASA contract NAS5-26555, and the NASA-CNES-CSA Far Ultraviolet Spectroscopic Explorer, 
which is operated for NASA by the Johns Hopkins University under NASA contract NAS5-32985. 

^Visiting at the Johns Hopkins University, Henry A. Rowland Department of Physics and Astronomy, 
Baltimore, MD 21218 
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absorption feature in the lUE LWP spectra. From the best-fit model solutions to 
the combined HST + FUSE spectral energy distribution, we find that the white 
dwarf is reaching a temperature T^ff ~ 45 — 55, 000 K in quiescence, assuming 
log{g) = 8.3 with a solar composition accreted atmosphere. The exact temper- 
ature of the WD depends on the reddening assumed and on the inclusion of a 
quiescent low mass accretion rate accretion disk. Accretion disk models alone fit 
badly in the FUSE range while, and if we take the distance to be a free param- 
eter, the only accretion disk model which fits well is for a discordant distance 
of at least several hundred pc and an accretion rate (~ 10^^Mq/jt) which is 
unacceptably high for a dwarf nova in quiescence. We discuss the implications 
of the white dwarf's temperature on the time-averaged accretion rate and long 
term compressional heating models. 

Subject headings: accretion, accretion disks - novae, cataclysmic variables - white 
dwarfs 



1. Introduction 

SS Cygni is a prototype dwarf nova, a subclass of the cataclysmic variable group of 
binary stars. Cataclysmic variables contain a white dwarf and a larger radius, less massive, 
donor star. The donor star fills its Roche lobe and its matter has gradually been cannibalized 
by the white dwarf. In dwarf novae the matter feeds into a disk around the white dwarf and 
continues to build up the disk mass until the disk reaches a critical temperature at which 
time the mass collapses onto the white dwarf and releases gravitational potential energy as 
radiation. This results in a roughly periodic increase in brightness of the system. 

SS Cyg is one of the best-studied dwarf novae. It has an orbital period of 6.6 hr 
(well above the period gap) and a ~ 50 da y recurrence time between dwarf nova outbursts 



(ICannizzo &: Mattel 1119921 : ICannizzo Ill993l ). These factors and its well-studied history have 
made it a key target of our ongoing studies to understand how disk accretion affects the 
white dwarf. Observations must be taken when the system is in quiescence because in such a 
state, the disk is least luminous which offers a favorable opportunity to detect the radiation 
of the white dwarf photosphere. However, the source of far ultraviolet light during the 
quiescence of dwarf novae remains controversial. This is especially true in dwarf novae with 
orbital periods above the CV period gap, which tend to have higher mass transfer rates, 
larger accretion disks and more massive secondary stars. Therefore, during quiescence it is 
possible that some portions of the remaining disk may still be optically thick and dominate 
the far ultraviolet (FUV) flux of the system. 
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The inclination of its orbit i s ~ 51° (IBitner et al. 1 120071 ) while its white dwarf mass 
has been reported to be 1.2Mq (IShafter I Il983l ) and secondary mass 0.7Mq. The system 
has an accurate paral lax measured with the Hubble Space Telesc ope fine guidance sensor 



(JHarrison et al. 



20001) yielding a parallax distance of 166l]^2 P^- 



Holm fc Polidan I fllQSSi ) 



first noted the possible detection of the underlying accreting white dwarf from the occasional 
appearence of Lya absorption in some lUE spectra obtained during quiescence. However, 
despite a number of attempts to measure the white dwarf's temperature and quantify its 
contribution to the FUV flux, the uncertain contribution of the accretion disk and other 
possible unidentified source (s) of FUV emission, the "second component" , rend ered any 
such measurements suspect ( JHolm Ill988t iLesniak fc Sion II2003I : iLong et al. Il2005[ ). 



Recently however, we retrieved archival FUSE and HST spectra, both obtained during 
different quiescent intervals of SS Cygni but with closely matching flux levels in the wave- 
length range where HST/GHRS and Ft/^i? overlap (~ 1150 — 1190 A ) . Also a new mass de- 
termination for the white dwarf has now been derived by lBitner et al. I (120071 ): O.81M0±O.18, 
significantly lower than the earlier value of 1.2M0. Using this lower mass and the parallax 
distance of 166 pc, we have reexamined the nature of the hot component with synthetic spec- 
tral models of high gravity photospheres and optically thick accretion disks. This assessment 
is presented below. 

In the next section we give details of the archival spectra, in section 3 we present our 
spectal modeling tools and method, results are presented in section 4 and discussed in the 
concluding section. 



Archival Observations 



The archival HST GHRS spectrum was taken on JD 2450352 (1996 September 26), 
in mid- quiescence, 17 days after SS Cygni reached quiescence and 21 days after the pre- 
vious outburst. The HST spectrum is a combination of two individual spectral segments 
Z3DV0204T (1150 A- 1435 A) and Z3DV0205T (1377 A- 1663 A) taken in ACCUM mode 
with the GHRS spectrograph with the G140L grating and LSA aperture (of size 1.74"). The 
two spectral segments, one at each grating setting, had an exposure time of 2176s each, and 
there was a 32.25 second time gap in an otherwise nominal exposure. The GHRS observa- 
tions were calibrated using the standard pipeline CALHRS, and were retrieved as VO-Tables 
using VOSpec. 

The archival FUSE spectrum was taken on JD 2452156.5 (2001 September 4), 18 days 
after the previous outburst, and 13 days after the system first entered quiescence. The 
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F^/S'i? spectrum (P2420101), was a combined spectrum of 8 separate exposures taken through 
the 30"x30" LWRS Large Square Aperture in TIME Tag mode. The total good exposure 
time was ~ 18ks, varying shghtly for each spectral channel. The observations were carried 
out mostly dur ing NIGHT time. The FUSE observations were processed with CalFUSE 
version v3.2.3 ( JDixon et al. 1120071 ). We follow the same procedure we used previously for 
the postproce ssing (co-addit i on, al ignment and weight of the spectral channels) of the FUSE 
data (see e.g. iGodon et al. I (120061 )). 



The combined FUSE + GHRS spectrum of SS Cyg is presented in Fig.l (see section 4 
for the model fit). The spectrum is characterized by strong emission lines originating in an 
optically thin regio n in the disk or above it, as inferred from their rotational broadening (see 



Long et al. I (120051 ) for a rigourous analysis of these lines in the GHRS spectra). In the very 
short wavelengths of FUSE the broad emission lines from N IV, S VI and H I (Ly5 and higher) 
merge together and produced an apparent rise of flux (< 950 A), which we do not attempt 
to model. In the fitting (section 4) we mask all the emission lines and these appear in blue 
in Fig.l. It is likely that around 1060-1080 A and 1110-1120 A some additional emission is 
present. 



3. Synthetic Spectral Modeling 



W e adopted model accretion disks from the optically thick disk model grid of lWade fc Hubeny 
(119981 ). In these accretion disk models, the innermost disk radius, Rm, is fixed at a fractional 
white dwarf radius of x = Rin/Rwd = 1-05. The outermost disk radius. Rout, "was chosen so 
that Tf,ff{Rout) is near 10,000 K since disk annuli beyond this point, which are cooler zones 
with larger radii, would provide only a very small contribution to the mid and far UV disk 
flux, particularly the FUV bandpass (~ 900 — 1700 A). The mass transfer rate is assumed 
to be the same for all radii. 

Theor etical, high gr avity, photospheric spectra were computed by first usi ng the code 
TLUSTY (JHubeny Ill988l ) to calculate the atmospheric structure and SYNSPEC (jHubeny fc Lanz 
19951 ) to construct synthetic spectra. We compiled a library of photospheric spectra covering 
the temperature range from 15,000 K to 70,000 K in increments of 1000 K, and a surface 
gravity range, \og{g) = 7.0 — 9.0, in increments of 0.2 in \og{g). 

The reddening of the system was taken from estimates list ed in the li t eratu re, determined 



from t he strength of the 2200 A interstellar absorption feature. IVerbunt ( 1987 ) and La Dous 



(Il99ll ) both give E(B-V)=0.04 while the more recent work of JBruch &: Engel I (Il994j ) gives 
E(B-V)=0.07. Both values are much smaller than the galactic reddening in the direction 
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of SS Cyg which is pretty large (~ 0.5) in agreement with the fact that SS Cyg is rather 
nearby with a distance of only 166 pc. The combined spectrum was de-reddened with the 
lUERDAF IDL routine UNRED assuming both E(B-V)=0.04 and E(B-V)=0.07. 

After masking emission hues in the spectra, we determined separately for each spectrum, 
the best-fitting white dwarf-only models and the best-fitting disk-only models using a x^ 
minimization routine. A x^ value and a scale factor were computed for each model fit. The 
scale factor, S, normalized to a kiloparsec and solar radius, can be related to the white dwarf 
radius R through: 



'A(ofes) 



SH 



\(model) ) 



where S = AnR d 



and d is the distance to the source. For the white d warf radii, we use the mass-radius relation 
from the evolutionary model grid of IWood I ( 119951 ) for C-0 cores. We combined white dwarf 
models and accretion disk models using a x^ minimization routine called DISKFIT. Using 
this method the best-fitting composite white dwarf plus disk model is determined based 
upon the minimum x^ value achieved, visual inspection of the model, consistence with the 
continuum slope and Lya region, and consistence of the scale factor-derived distance with 
the adopted trigonometric parallax distance. 



4. Synthetic Spectral Fitting Results 

We used an accretion disk model, a white dwarf photosphere, and a combination of 
both in our analysis. The disk models used were optically thick, steady state accretion 
models with solar abundances. They are considered a reasonable first approximation to the 
spectral shape of the disk in quiescence. In the disk models, we first adopted an inclination 
of 41° and 60° directly from the grid of models of IWade fc Hubeny I (119981 ). Only after a 
satisfactory best fit was obtained with the correct (parallax) distance, did we then generate 
disk models (using TLUSTY, SYNSPEC and DI SKSYN) with an inclination of 50°, closer 
to its derived value of 51° ±5 (JBitner et al. 1120071 ). The photosphere models were generated 
using TLUSTY and SYNSPEC. In the following we used the parallax distance of 166 pc, 
and the mass of 0.8Mq (log(f7) = 8.3) (JBitner et al. 1120071 ). We assumed reddening values 
of E(B-V) = 0.04 and 0.07. A summary of the model fits is given in Table 1. 

We started our fitting with single accretion disk models alone assuming both E(B- 
V)=0.04 & 0.07. We ran disk model fits and found that the best fit (lowest x^) had a mass 
accretion rate far too large for quiescence and a distance much larger than 166 pc. 

Next, we tried single white dwarf atmosphere models, first dereddening assuming E(B- 
V)=0.04. Since almost all the lines are in emission (possibly from an optically thin region 
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in the disk and/or corona), it makes it difficult to assess the rotational velocity broadening 
and chemical abundances based on the absorption lines. Nevertheless, for all WD models 
we assumed solar abundances and a canonical projected rotational velocity of 200 km/s, 
and checked that the value of the x^ did not depend on Vrot sin ^ as long as it was a few 
hundred km/s. The best fit least x^ WD model has a temperature Te// = 40, 000 K and 
X^ = 1.637, but with a distance of only 138 pc. A model with Tg// = 47,000 K gave the 
right distance with a slightly larger x^, namely 1.990. Next we ran single WD model fits 
assuming E(B-V)=0.07. The x^ we obtained increased slightly over the E(B-V)=0.04 best 
fit models. 

Last, we explored whether the fitting could be improved if we combined a white dwarf 
model with an accretion disk model. We found that some of the white dwarf plus disk 
combinations yielded distances close to 166 pc with a lower x"^- We have summarized some 
of these combination fits in Table 1. For E(B-V)=0.04, the best WD+disk fit leading to a 
distance in agreement with observed parallax is for a WD with T = 46, 000 K, a disk with 
a mass accretion rate M = 1 x lO^^°M0/yr, i = 50°, where the WD contributes 88% of the 
fiux and the disk contributes the remaining 12%. This model fit is presented in Figure 1. 
We then carried out the same fitting but this time assuming E(B-V)=0.07 and found similar 
results for the disk but with a higher WD temperature: T^ff = 55, 000 K. 

Therefore, we conclude that the dominant source of the far UV radiation between 912 
A and 1660 A is an accretion-heated white dwarf photosphere with T^ff si 45 — 55, 000 K, 
\og{g) = 8.3 (depending on the reddening value). 



5. Conclusions 

We have presented evidence from our model fitting analysis of the combined FUSE + 
HST/GHRS spectrum of SS Cyg during quiescence that the source of the far ultraviolet 
continuum and absorption line radiation is the white dwarf's photosphere. The disk models 
that best fit the spectral data yield unreasonably large distances, at odds with the HST 
FGS parallax and indicate accretion rates far too high to be associated with dwarf nova 
quiescence. The photosphere models give effective temperatures of 45,000 to 55,000 K for a 
reddening of 0.04 and 0.07 with the inclusion of a low mass accretion rate disk in agreement 
with the quiescent state. Unfortunately, the lack of a clear detection of absorption lines due 
to accreted metals and helium in the WD atmosphere prevents a determination of both the 
rotational velocity of the white dwarf and the abundance of metals in its accreted atmosphere. 

Our derived temperature for the white dwarf in SS Cygni is well above the presently 
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estimated average temperature (~30,000 K) for white dwarfs in dwarf novae above the CV 
period gap. Compared with the temperatures of white dwarfs in dwarf novae whose orbital 
peri ods are close to the period of SS Cygni, the white dwarf in TT Crt is cooler (40,000 



K; (ISion et al. II2008I )) while the white dwarf in Z Cam is hotter (57,000 K; ( iHartley et al. 



20051 )). If the accretion rate scales with the orbital period, then the temperatures should be 
comparable. It is interesting that the white dwarfs in dwarf novae with Porb < 360 mins are 
cooler than 40,000 K while the hottest white dwarfs in dwarf novae are found at Port > 360 
mins. 

At Porb = 6.6 h, the Tg// of the WD in SS Cygni lies within the range expected from 
compressional heating for an average M , < M >, obtained from typic al interrupted magnetic 
braking laws for white dwarf masses between O.6M0 and I.OMq (JTownsley &: Gansicke 



2OO9I ) . A linear extrapolation to Porb, of the predicted T^fj for Porb = 6.6 h, corresponds to 



an average mass transfer rate of < M >~ 10 ^Mq/jt which is at the high end of the range of 
mass transfer rate s associa.ted w ith the nova-like variables, as deter mined from their optica l 



disk luminosities ( Warner 



I995I ). Interestingly, the braking laws of lAndronov et al. I (120031 ) 



and llvanova fc Taam I ( 20041 ) either fall drastically short or exhibit a downturn, respectively, 
of yielding the < M > implied by the WD T^jf, while the iHowell et al. I (1200 ll ) law steeply 
increases at constant Porb ~ 5 h. One cannot rule that other significant sources of heating 
of the white dwarf besides compression are operating such as possible nuclear burning. It 
seems clear that more white dwarf temperatures are needed in dwarf novae and nova-like 
variables at long Porb before definitive conclusions can be drawn. 
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Table 1. Synthetic Spectral Model Fits 
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Fig. 4. — The combined spectrum of SS Cyg (in red) dereddened assuming E(B-V)=0.04 
is shown with a composite WD+disk synthetic spectrum (in sohd black). The WD model 
(dotted line) has My^a = O.SMq (corresponding to Log{g) = 8.3), T^ff = 46, 000 K, solar com- 
position and a projected rotational velocity of 200km/s. The dashed line is the contribution 
from the accretion disk model. The disk model has a mass accretion rate M = 10~^°MQ/yr 
and an inclination i = 50°. The x^ obtained is 1.258 and the resulting distance is 173 pc. 
The WD contributes 88% of the FUV flux, while the disk contributes the remaining 12%. 
The excess of flux around 1060-1080 A , and 1110-1130 A might be due to some emission 
from Siv (1063, ~1073 A), Siiv (1066A), Sim (~1108-1113 A), and Si IV (1122.5, 1128.3 



